Charmed Baryonium 



G Cotugno t , R Faccinit.K, AD Polosa^, C Sabelli t ' t 
1 INFN Roma, Piazzale A. Moro 2, Roma, 1-00185, Italy 
Department of Physics, Universita di Roma, 'La Sapienza', Piazzale A. Moro 2, Roma, 1-00185, ttaly 

We re-analyze the published data on the Y(4630) -¥ A C A C and the F(4660) -¥ ip(2S)-Kn 
with a consistent Ansatz and we find that the two observations are likely to be due 
to the same state Y B with My b = 4660.7 ± 8.7 MeV and Ty B = 61 ± 23 MeV. Un- 
der this hypothesis and reanalizing also the e + e~ — > J/ipinvj ISR spectrum we extract 
B(Y B -»■ A C A C )/B(Y B -> V(2S>tt) = 25 ± 7, B(Y B -> J/ipirK)/B(Y B V(2S>7r) < 0.46 
@ 90% C.L., B(Y(4350) ->■ J/-07r7r)/S(F(435O) -> ip(2S)iriv) < 3.4 x 10~ 3 @ 90% C.L. and 
B(Y B — > iP(2S)ct)/B(Yb — > ip{2S)fo) = 2.0 ±0.3. These conclusions strongly support the hypothesis 
of Y B being the first observation of a charmed baryonium constituted by four quarks. From the 
analysis of the mass spectrum we show that F(4350) and Y B are respectively consistent with the 
ground state and first radial excitation of the I = 1 state. 

PACS: 12.39.-x, 12.39.Mk, 13.75.-n 



Introduction. The search of the so called baryonia, exotic states originally thought to appear in nucleon- 
antinucleon systems, has a rather long and sometimes controversial history. Early work in this field can be 
found in [IJ and quite recently there have been new interesting experimental indications as discussed e.g. in [2j. 
In this letter we focus on higher mass scales, namely we refer to the hidden charm sector. We identify as a 
charmed baryonium a narrow structure, which we call Yb(J pc = 1 ), with a decay pattern dominated by the 
A C A C baryon-antibaryon mode and compare a set of other observablcs with the expectations of a tetraquark 
model. Many different interpretations have been proposed in previous works on the nature of the 1 neutral 
states in the 4 GeV region. The conventional charmonium assignation has been discussed in Q. Besides the 
charmonium hybrid interpretation Q and the tetraquark exlpanation [Bj], the possibility that K(4260) could be 
a threshold effect [6] or a charm mesons molecule [7[ have been considered too. Recently [8( and [9( proposed 
the F(4660) to be a D*D\ and a ip fo bound state respectively. 

Yb from Belle data. The existence of Yb stems from our re-analysis of data from the Belle experiment 
on the decay of the y(4660) resonance in i/j(2S)Tnr [lfj and of the y(4630) in A C A C [llj. We fit the invariant 
mass spectra in Fig. [T] with a binned likelihood, adopting a consistent signal model: a relativistic Breit-Wigner 
with comoving width, as detailed in Ref. (l2l |. Background is parameterized with a second order polynomial 
multiplied by the phase space. 

The individual fits to the charmonium and baryonic modes return My B = 4661 ± 9 MeV, ry B = 61 ± 23 MeV 
with x 2 /d.o.f. = 7/20 and M Yb = 4661 ± 14 MeV, I> B = 63 ± 23 MeV with x 2 /d-0.f. = 51/35, respectively. 
The two results are consistent, strongly supporting the hypothesis that the two structures are evidences of the 
same resonance. In the following we will use the results of the first and more accurate, fit as reference. From 
the same fits we also extract 

B(Y B -> A C A C ) 

= 25 ±7, (1) 



B(Y B -> V(25*)7r+vr- 

a result which highlights a strong affinity of the Yb to the baryon-antibaryon decay mode. It is to be noted 
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FIG. 1: Fit of the cross section from the data of Belle Collaboration [TO, HH for the process e + e~ — > ip{2S)ir + n" (left) 
and e + e~ —¥ A^A~ (center) and from the BaBar Collaboration [lj| for the process e + e~ —¥ J/i/j-k + tt~ (right) . The 
black line represents the fit results, the red one shows the polynomial background. 



1 



that omitting interference, the ratio of the peak cross sections is about 1 1 . 

Another puzzling aspect of this meson is the absence of the decay mode into J/?/>7r + 7r~. In order to try to 
interpret this observation, we also extract, with the same signal model, an upper limit on the data from BaBar on 
this decay (see Fig. [T][l3j]) and extract, with a bayesian approach, B(Yb — > J/^tttt)/B(Yb — > ^(2S)tttt) < 0.46 
@ 90% C.L. The presence of the nearby Y(4260) was properly accounted for. 

By fitting consistently the i/j(2S)ir + Tr~ spectrum in Fig. [1] we also extract the parameters of the Y(4350) 
resonance. We obtain M = 4353 ± 13 MeV and T = 118 ± 25 MeV. We also find that £(Y(4350) -4 
J/V>7T7r)/23(T(4350) -> i/i(2S)mr) < 3.4 x 10~ 3 @ 90% C.L. 

Finally we considered the 7T7t invariant mass spectrum as measured by Belle in the mass region of the Yb [13] 
(see Fig. [2]). In absence of background (see Fig. Q}, we fit for the sum of two relativistic Breit-Wigner functions 
describing the a and /o(980). The resulting fitted masses and widths are m a = (714 ± 77) MeV, T a = 
(499 ± 176) MeV and m fo = (955.3 ± 4.5) MeV, r /o = (14 ± 15) MeV, in agreement with the standard 
knowledge about the light scalar mesons. Since these mesons are the best currently known candidates for 
tetraquarks [l4| , the four-quark nature of the Yg is naturally suggested. Furthermore we can extract the ratio 

B(Y B ip{2S)a{600))/B{Y B -> V(2S)/ (980)) = 2.0 ± 0.3. (2) 
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FIG. 2: Fit of the dipionic mass distribution in the mass region for ip(2S)'K + ir~ : 4.0 GeV < m^+^-^nm < 4.5 GeV [llj 
including the contributions of a and /o(980). Removing the a, the lower part of the spectrum gives a sensibly worse fit. 

Interpretation of Yb- The dominance of decay modes with two baryons or with scalar mesons are a 
strong indication of a [cq] [eg] diquark-antidiquark structure of the Yb meson. The diagram A in Fig. [3] is a 
representation of a [eg] [eg] tetraquark where the two diquarks are connected by a string neutralizing their color. 
Both diquarks have positive parity in the 'single mode configuration', where the constituent quarks are unexcited 
with respect to one another. A J PC = state must therefore have odd values of I. The 1 = 1 centrifugal 
potential favors configurations with the diquarks at an higher average relative distance with respect to the i = 
case, stretching the color string which eventually breaks as in diagram B of Fig. [3] (baryon-antibaryon decay). 
Breaking the color string in two points rather than in one, as shown in C or in D in Fig. [3] allows either the 
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FIG. 3: Phase space allows A -»■ B (Y -> A C A C ), A — > C (Y — > DD*), A -> D (Y -> i/'(2S , )/o(980)) 
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DD* or the ipfo decay. Actually also I = 3 would be possible assuming the total spin of the^diquarks to be S = 2 
with S = S„ + S* . Although an analysis of the mass spectrum made along the lines of [l7| would exclude this 
possibility, we will show that a different implementation of the orbital angular momentum dependency of the 
Hamiltonian would still allow it. We aim at estimating the mass of the Ys using a non-relativistic constituent 
quark model combined with a hadron string model. Using the fitted masses of constituent diquarks we will find 
the spectrum as a function of the orbital angular momentum and radial excitation quantum numbers. 

As for the isospin assignment of the tetraquark wavefunction the observed decays seem to favor I — 0. On 
the other hand at the mass scale of the charm quark we could assume that light quark annihilation diagrams are 



suppressed so that the mass eigenvectors align to the quark mass basis |15|. In the latter case one could have 
isospin breaking effects as is the case for the A(3872) which turns out to decay in J/tp p and J/ipuj at almost 
the same rate. In a tetraquark interpretation one could explain such decay dynamics making the hypothesis 
that there are two states [16| of the kind X u = [cu] [cu] and = [cd] [cd] with a difference in mass of the order 
of 7 x cos(26») MeV, 6 being the mixing between X u and X d ~ 20° in the case of X(3872)). From available 
data we still do not know if this is the case for the Yg particle although this is not affecting the following 
considerations. 

Angular excitations. The tetraquark picture with two diquarks connected by a color string allows to 
estimate the mass of the Yb and its orbital angular momentum excitation level. Considering the generalization of 
the Chew-Frautschi formula discussed in Selem and Wilczek [1S|] we find in the mass limit of Mu) (M = mi = 7712) 
very large with respect to a/2ir 1 

f * 2M +(rdb^ (rf)2/3 <3) 

where £ is the energy of a system made up of two diquarks of mass M = held together by a string 

characterized by a tension a and spinning with some angular momentum I. The higher t the better this model 
works as shown in Regge trajectories phenomenology. Taking a, fitted from Regge slopes of non-charmed 
hadrons, to be a ~ 1 GeV 2 and M = 1933 MeV . we have £ w 4311 MeV for I = 1 and £ w 4793 MeV for 
I = 3. This means that if spin-spin interactions are neglected, which is a good approximation at £ ^ 0, and 
if we consider no radial excitation, the I = 3 assignation might be favored while the t = 1 state would better 
match the Y(4350). 

To check the effectiveness of Eq. ((3]) in estimating the mass of known charmonium resonances having different 
I it is worth noting that, for M — 1500 MeV, we get 3500 MeV for the average mass value of the Xcj system 
(experimentally located at ~ 3493 MeV), which has I = 1, and 3770 MeV for the I = 2 case consistent with the 
D-wave nature of the ^(3770) [2lj. The fine structure of the XcJ system, with levels shifted from the average 
mass value, is due to spin-orbit and tensor interactions. 

Radial excitations. Let us assume that Yg is a diquark-antidiquark bound state in a Charmoinium-like 
potential. The difference in mass between the n r = and the first radial excitation n r = 1 can be extracted 
from the available data on Xbj{2P), Xbj(^P) states. The mass difference AM = M Xb r 2 p\ — M Xh n P ^ ~ 360 MeV 
is expected to be roughly the same for charmed states, where instead no data are available for £ = 1 states. 
Therefore, neglecting spin-orbit terms and spin-spin interactions, the (n r = 1,^ = 1) state would be located at 
4671 MeV, remarkably close to Yg(4660), whereas the (n r = 0,£ = 1) state would be located at 4311 MeV, 
close to the Y(4350). We can make an estimate of the mass spectrum also using the constituent quark model 
discussed in [La |. If we neglect spin-spin and spin-orbit splittings and assume I = 1 we find My = 4361 MeV if 
no radial excitation is taken into account, and My B = 4721 MeV if n r = 1. Therefore we might conclude that 
Y B (4660) = Y(n r = 1,£=1) and Y(4350) = Y(n r = 0,£= 1). 

Decays. To refine the model we investigate whether the decay patterns of the Y(4350) and Yg(4660) are 
consistent with the assigned quantum numbers: from data we see that both Y(4350) and Yg(4660) prefer to 
decay in tJj(2S)tttt rather than in J/ipirir even if the latter mode is phase space enhanced. The measured spectra 



1 We perform an infinite mass limit of the the equations of a rotating relativistic string with tension T = a/2-n and with two 
masses m\ , m,2 attached at its ends, namely 19j 



mi , ni2 , a f ur i dv a dv 



r ur i dv a r 

yjl - (turi) 2 ' y/1 — (W 2 ) 2 ' 2-n-oj J •/% - v' 2 2-klo J y/l — v 2 
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show that the tttt pair comes from phase space (or, equivalently a a) apart from the case of the Yb (4660) where 
there seems to be a 30% component due to /o(980) (see Eq. 2). We therefore attempt an explanation of the 
observed Y's decay pattern describing the S'-wave transition (ip(lS,2S)s\Y) as 

(i()(ri,q)a(k)\Y(e,p)} =ge-r} (4) 

where a = <r, fo- We find the eigenfimctions ^ of the linear part of the Cornell potential V(r) = —k/r + r/a 2 
(k = 0.52, a — 2.3 GeV -1 ) for the ip(lS, 2S) and the charmonium-like Y, namely ^ = Y^ m R nr i. Interpreting 
the Y's as charmonium-like bound states made up by a diquark and an antidiquark we can write 

gcx J d 3 r rt w (r)R (Y) (r) (5) 

We are interested in estimating the following ratio of decay widths 

r a (Y) = T(Y -> ^(lS){mr) a )/T{Y -> ^(25)(7T7r) a ) (6) 

We find the values reported in the following Table (the square of the matrix element in Eq. §5§ weights the 
three-body phase space where the two pions are the decay products of an intermediate 2 scalar resonance; we 
include only those quantum number combinations which provide mass values close to the experimental ones) 
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We report separately the a and /o contributions. The results show that the radial excitation would explain the 
preference of Yg tetraquark to decay into ip(2S)a whereas with this simple model we cannot explain why the 
Y(4350) seems to prefer the ip(2S) channel. This aspect remains puzzling. 

Conclusions . We reanalized the published data on the Y(4630), Y(4660) and Y(4350) and obtained that 
the first two resonances correspond to a single state, the Yb (4660). The large baryon decay rate and the decay 
into t/>(25)/o(980) make the Yb an excellent candidate for a [cd][cd] diquark-antidiquark bound state suggesting 
that Yb represents the first observation of a what we might call a charmed baryonium. The fact that the 
mass of the Yb is sensibly higher than the baryon-antibaryon threshold excludes any baryon molecule picture. 
Assuming that the interaction between the two diquarks has the same potential as that between quarks in the 
charmonium we estimate the mass spectra for the ground state and the first radial excitation for two values of 
I and find that Yb best fits the first radial excitation of I = 1. With the same assumptions we show that the 
Y(4350) fits the corresponding radial ground state and we show that the absence of the decay into J/ipmr of 
the Yb (4660) seems compatible with its tetraquark interpretation. 
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